Poisson commutator-anticommutator brackets for ray tracing 
and longitudinal imaging via geometric algebra 
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Abstract. We use the vector wedge product in geomet- 
ric algebra to show that Poisson commutator brackets 
measure preservation of phase space areas. We also use 
the vector dot product to define the Poisson anticom- 
mutator bracket that measures the preservation of phase 
space angles. We apply these brackets to the paraxial 
meridional complex height-angle ray vectors that trans- 
form via a 2 x 2 matrix, and we show that this transforma- 
tion preserves areas but not angles in phase space. The 
Poisson brackets here are expressed in terms of the coeffi- 
cients of the ABCD matrix. We also apply these brackets 
to the distance-height ray vectors measured from the in- 
put and output side of the optical system. We show that 
these vectors obey a partial Moebius transformation, and 
that this transformation preserves neither areas nor an- 
gles. The Poisson brackets here are expressed in terms 
of the transverse and longitudinal magnifications. 



1 Introduction 

a. Poisson Brackets. In two dimensional space, the 
Poisson bracket p] is defined by the commutator or the 
Jacobian 
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(1) 



where we used the subscript notation of Hand and 
Finch 2 . If this bracket is unity, then the mapping from 
the object parallelogram dq A dp to its image dF A dG is 
a symplectic (area-preserving) map. 3 

So we ask: if the Poisson bracket commutator in 
Eq. P is a measure of area preservation, can we also 



define a similar bracket for the preservation of the inte- 
rior angles of the differential parallelograms dq A dp and 
dF A dG? Since this problem appears not to be taken 
before, we shall propose the following anticommutator 
bracket measure: 



_ d£d£ dGdG 

q ' p dq dp dp dq ' 



(2) 



which we shall show later to be based on the definition 
of the dot or inner vector product. If dq _L dp and the 
bracket in Eq. ^ is zero, then dF J_ dG. 

In geometric algebra, the wedge (outer) and dot (in- 
ner) products of two vectors a and b are related by their 
geometric (juxtaposition) product [1]: 



ab 



b + aAb. 



(3) 



The dot product is the same as in the standard vector 
analysis; the wedge product is the same as the Grassman 
outer product [5 . Doran and Lasenby used the wedge 
product of two vectors to express the Poisson commuta- 
tor bracket commutator in Eq. {1} as 6J 



[F,G] = (Vf AVG) • J, 



(4) 



where J is a bivector (two form or oriented area); how- 
ever, they have not used the dot product to the define 
the bracket's anticommutator counterpart. 

b. Paraxial Optics. In paraxial meridional ray op- 
tics, a light ray is described by the height angle vector 
(x, not), where x is the height of the light particle from 
the optical axis, n is the refractive index of the medium, 
and a is the inclination angle of the direction of propaga- 
tion of light. The input (unprimed) and output (primed) 
light rays are related by the ABCD system matrix: [7] 
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Though this linear transformation is known to preserve 
the area dx A d(na) , we shall show that this transforma- 
tion does generally preserve the angle between dx and 
d(na). 

For imaging systems, the position of the object is given 
by (S,x) and that of the image by (S',x'), where S — 
s/n is the reduced distance of the object measured from 
the input side of the optical system, and S' = s'/n' is 
that of the image measured from the output side. These 
parameters are related by the partial Moebius relations 

M 11 S + M 2l 

M 12 S - M 22 ' W 



M 12 S-M 22 ' y ' 

where M\x, Mi 2 , M 2 i, and M 22 are constants. We shall 
show that the transformation described by Eqs. ([6]) and 
([7|) is not symplectic, because the area dS A dx is not 
preserved. We shall also show that this transformation 
does not preserve angles. 

c. Outline. We shall divide the paper into five sec- 
tions. The first section is Introduction. In the second 
section, we shall present geometric algebra, we we shall 
discuss products of vectors and of complex vectors. In 
the third section, we shall summarize the matrix equa- 
tions in paraxial meridional ray tracing in complex vec- 
tor formalism we proposed in our previous paper [8]. We 
shall discuss the linear height-angle transformation and 
the Moebius-like distance-height transformation. In the 
fourth section, we shall compute the Poisson commuta- 
tor and anticommutator brackets of these transforma- 
tions. We shall show that the height-angle transforma- 
tion preserves areas but not angles in phase space, while 
the distance-height transformation neither preserve areas 
nor angles of longitudinal objects. 

2 Geometric Algebra 

The Clifford (geometric) algebra Cl 3 fi is a group[9] alge- 
bra over the field of real numbers. The generators of the 
group are the three vectors ei, e 2l and e3 that satisfy 
the orthonormality relation [10J 

ejefe + e k e 3 = 25 jk , (8) 

for j,k = 1,2,3. That is, 

ef = 1, (9) 
<w = <V <v (10) 

for j ^ k. Equation ([9]) is called the normality axiom: 
a vector with a unit square has a unit norm or length. 
Equation (Tl0|) is the orthogonality axiom: perpendicular 
vectors anticommute. 



2.1 Vector Products 

Let a and b be two vectors spanned by the three unit 
spatial vectors in 01%$: 

a = aiei + a 2 e 2 + a 3 e 3 , (11) 
a = aiei + a 2 e 2 + a 3 e 3 . (12) 

By the orthonormality axiom in Eq. l[8]). we can show 
that the geometric (juxtaposed) product of these two vec- 
tors is 

ab = ab + aAb, (13) 

where 

ab = aih + a 2 b 2 + a 3 b 3 , (14) 
aAb = eie 2 (aib 2 ~ a 2 bi) + e 2 e 3 (a 2 b 3 - a 3 b 2 ) 

+ e 3 ei(a 3 6i - aib 3 ) (15) 

are the inner (dot) and outer (wedge) products of the 
two vectors [Tl]. 

To geometrically interpret the dot and wedge products 
of two vectors, let us define a as a vector along ei and 
b as the vector along ei rotated counterclockwise about 



the vector e 3 by an angle tf>: 

a = ae%, (16) 
b = ei&cos (p + e 2 bsin<f>. (17) 

The product of a and b is 

ab = ab cos cf> + eie 2 ab sin cf>. (18) 

Separating the scalar and bivector parts, we get 

a-b = abcoscf), (19) 
aAb = ie 3 ab sin cf> = eie 2 absmcj). (20) 



Thus, the scalar a • b is the product the magnitude a 
of vector a and the component of b cos <p of the vector b 
along a; the bivector a A b is proportional to the area of 
the parallelogram defined by vectors a and b. (Fig. [T]) 




a 



Figure 1: An oriented area aAb defined by the 
vectors a and b. The angle between the two vec- 
tors is cf>. 
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The wedge product aAb in Eq. (|15|) may be expressed 
in determinant form as 



Ol 


6l 






l>2 




a 3 b 3 




b- 2 


eie 2 + 




h 


e 2 e 3 + 


Ol 61 


«2 




a 3 



a Ab = " eie-2+ , o-jiv, " e3ei. 

Using the properties of determinants, we can show that 



aAb = b A a, 
aAa = 0. 



(22) 
(23) 



The first equation states that the orientation of the di- 
rected area defined by two vectors flips if the vector fac- 
tors are interchanged. The second equation states that 
no area can be defined by a vector and itself. 

The dot product is a grade-lowering operation; the 
wedge product is a grade-raising operation. That is, 



a b = (ab)|i_i| = (ab) , 
a A b = (ab)| 1+1 | = (ab} 2 , 



(24) 
(25) 



where the angle brackets (ab}u denotes the extraction 
the g— vector part of ab. In other words, the grade of 
a • b is the sum of the grades of a and b; the grade of 
a A b is the difference of the grades of a and b. 



where 



a x b 



(a 2 b 3 - a 3 6 2 )ei + (a 3 6i 
+ (ai& 2 - a 2 6i)e 3 . 



ai6 3 )e 2 



(32) 



Geometrically, the bivector a A b is the oriented plane 
with a x b as its associated normal vector. (Fig. [2]) 



a x b 
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Figure 2: The oriented plane aAb and its normal 
vector a x b. 



2.2 Hodge Map and Spatial Inversion 

Let us denote the product of the three orthonormal vec- 
tors ei, e 2 , and e 3 by z: [T2] 

i = eie 2 e 3 . (26) 

The geometric interpretation of i is an oriented volume. 

We chose the notation i because the trivector behaves 
like a true unit imaginary scalar. That is, its square is 
-1, 

i 2 = e 1 e 2 e 3 e 1 e 2 e 3 = e^e 2 e 3 e 2 e 3 

= -e?e|e§ = -l, (27) 

by the orthonormality axiom in Eq. (|8|) , and that it com- 
mutes with vectors, 

iei = e 2 e 3 = e\i, (28) 
ie 2 = e 3 e! = e 2 i, (29) 
ie 3 = eie 2 = e 3 i. (30) 

Notice that the trivector i transforms a vector into a 
bivector (and a bivector into a vector). This transforma- 
tion is called a Hodge map or a duality transformation. 

Using the relations in Eq. (|28| to ([30)1 . we can show 
that the wedge product aAb is related to the cross prod- 
uct a x b by a Hodge map: 

aAb = j(axb), (31) 



Substituting the Hodge map relation in Eq. (|31[) back 
to the juxtaposed product expansion in Eq. (|13p . we 
get[T3] 

ab = a- b + t(a x b). (33) 

This equation is similar to the Pauli identity in Quan- 
tum Mechanics, except for the absence of the Pauli 
a— matrices. 

In general, every element A in CZ 3j o may be expressed 
as a sum of a scalar, a vector, a bivector (imaginary vec- 
tor), and a trivector (imaginary scalar): [13] 

A = A Q + Ai+iA 2 + iA 3 . (34) 

We call such a sum a cliffor. The spatial inverse of the 
cliff or A is[T4"] 

i f = A - A x + iA 2 - iA 3 . (35) 

That is, the spatial inversion operator (t) flips the signs 
of vectors and trivectors, but leaves scalars and bi vec- 
tors unchanged. The other term for spatial inversion is 
automorphic grade involution [15], 

Other properties of the spatial inversion operator are 
as follows: 

(i + B) f = ti + B\ (36) 
(AB)t = (37) 
(it) f = A. (38) 
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These properties can be derived from the definition of the 
spatial inverse of A in terms of the unit vector eo = e4 
in CZ4 o that is orthogonal to ei, ea, and e3:[14"] 



ie = e i t . 



(39) 



Notice that spatial inversion does not flip the order of the 
products, unlike in the case of the reversion [16] (though 
they share the same dagger notation). 

2.3 Complex Vector Products 

Let r and f' be complex vectors: 



a + ib, 
c + id. 



The product of f and r ' is 

rr' — (a + ib)(c + id). 



(40) 
(41) 



(42) 



Separating the scalar, vector, bivector, and trivector 
parts of Eq. (|42|) , we get 



(ff')o — a ■ c — b • d, 
(ff')i = -axd-bxc. 
(ff') 2 = i(axc-bxd), 
i(a ■ d + b • c). 



(rr'h 



Take note of Eq. (|44|) . 

On the other hand, the product of f and r'> is 

ff'^ = (a + ib)(— c + id). 



(43) 
(44) 
(45) 
(46) 



(47) 



Separating the scalar, vector, bivector, and trivector 
parts of Eq. (|47]) , we get 



(48) 
(49) 
(50) 
(51) 



Take note of Eq. (|48|) . 

To clarify the geometric interpretations of Eqs. (|44|) 
and (14"5]) , let us go to the complex vector phase space 
and replace the imaginary parts by vectors. (Though this 
procedure is strictly not allowed, it is still pedagogically 
instructive.) That is, we write 



(rr'T)o = — a • c — b • d, 

(rr'')x = — a x d + b x c, 

(ff'^) 2 = i(— a x c — b x d) 

(ff'T) 3 = i(a-d — be). 



r = c 



b, 
d. 



The product of the vectors r and r' is 
rr' = (a + b)(c + d). 



(52) 
(53) 



(54) 



Separating the scalar and bivector parts of Eq. (|54p yields 

rr' = ac + ad + bc + b-d, (55) 
i(r x r') = i(a xc + axd + bxe + bx d).(56) 

Note that the scalar part is proportional to the cosine of 
the angle between r and r', while the magnitude of the 
bivector part is the area defined by r and r'. 
If we assume that a II c, b II d, and a J_ b. then 



ad = bc = axc = bxd = 0. 

Substituting these results back to Eqs. ([55)) and ([5 
obtain 

r • r' = a • c + b • d, 
i(r x r') = i(axd + bxc). 



Notice that r-r' is the negative of that (ff't) in Eq. 
and that the area |r x r'| is the same as the magnitude 
of (ff')i in Eq. (|4"4"|) . Later, we shall use these results 
later when we define the Poisson commutator and anti- 
commutator brackets in complex vector spaces. (Fig. [3]) 



(57) 
we 

(58) 
(59) 




Figure 3: The oriented area (a + b) A (c + d). 
This only serves as a way of visualizing the phase 
plane defined by the complex vectors f = a + ih 
and r — c + id. 



3 Matrix Optics 

3.1 Height- Angle Relations 

In matrix optics, a height-angle ray vector f may be de- 
scribed by a complex column vector, which we shall in- 
terpret in terms of the orthonormal vectors ei pointing 
upward and e 2 pointing out of the paper: [17] 



x 
ina 



xe\ + inae2, 



(60) 
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where x is the height of the light particle with respect to 
the optical axis e3 pointing to the right, n is the refractive 
index of the medium containing the light particle, and 
a is the counterclockwise paraxial angle of inclination of 
the direction of propagation of the light particle (Fig. [4]). 
Note that the imaginary number i = eie2e3 is a trivector, 
and ie.2 = e^ei is the zx— plane (bivector) containing the 
angle na. 

Let M be a complex system matrix, 



A -iC 
-iB D 

and let its determinant be unity, 

I M I = AD + BC = 1. 



(61) 



(62) 



Notice that because of the presence of the imaginary 
numbers, the determinant is not a difference, as given 
in the literature, but a sum. 



Separating the ei and ie 2 components of Eq. (f64|) . we 
get 



= Ax + Bna 
= —Cx + Dna. 



(65) 
(66) 



Except for the sign of C, these equations are similar to 
those in the literature. 

3.2 Distance-Height Relations 








x 2 




Xi 





tt2 



X3 



J 



Si 



*2 



Figure 4: A paraxial ray with height x\ and incli- 
nation angle ai moves to the right by a distance si 
until it hits an refracting surface. The height of the 
ray becomes X2 and its inclination angle changes 
to Ct2- 



When a light particle, characterized by its height-angle 
ray vector f, enters the (left) input side of an optical sys- 
tem described by the right-acting matrix M, the height- 
angle ray vector of the light particle changes into f at 



the (right) output side. In other words, 
f' = M T • f = f ■ M. 



(63) 



That is, Q3] 



fx' 
\ina 
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-iB\ 
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D ) 
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( A 
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ma ] 


\-iB 


D 



(64) 



Figure 5: In an imaging system, rays leaving a 
point source converge to a point. The distance of 
the object to the left input side of th black box is 
S, while that of the image from the right output 
side is S' . 



Let r be the position of an object point measured from 
the input side of the black box and let r' be the image 
point measured from the output side (Fig. [5]): 



where 



r = —Sea + xei, 
r' = 5"e 3 + x'ei, 



S = s/n, 
b = s /n 



(67) 
(68) 

(69) 
(70) 



are the reduced object and image distances, while x and 
x' are the object and image heights. At these two points 
the complex height-angle ray vectors are given by f and 
f , respectively. 

Let us decompose the system matrix M as the prod- 
uct of a propagation matrix S, a box matrix M|, ox , and 



where 



M.c 



M = M s M bo:E Ms< 



1 

-iS 1 



Mil 
-iM 2 i 

1 

-iS' 1 



-iM 12 
M 22 



(71) 

(72) 
(73) 
(74) 
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Notice that since M, Ms, and Mg/ have a unit determi- 
nant, then Mbox must also have a unit determinant: 



|M box | = M n M 22 + M 12 M 21 . 
The transpose of Eq. (fTTj) is 



1 T = M|,M^ ox M|. 



(75) 



(76) 



Taking the transpose of the matrix M in Eq. (ftJTj) and 
substituting the result in Eq. (|76p . we arrive at 

A = M n -M 12 S', (77) 

B = M 21 + M 22 S' + M 11 S-M 12 SS I , (78) 

C = Mia, (79) 

D = M 22 -M 12 S, (80) 

after carrying out the matrix products and separating 
the matrix coefficients. 

Now, if r and r' in Eqs. (f6"7) and are positions of 
the object and image points, then the object and image 
heights x and x' are constants, so that the output ray 
angle a' depends only on the input ray angle a: 



(81) 



a = a'(a). 



Thus, taking the partial derivative of Eqs. (f65|) and ([66 
with respect to na,we get 



= B, 



d(nW) 
d(na) 



= D. 



(82) 
(83) 



Our interest in this paper is only with the first equation. 

Substituting Eq. ^ into Eqs. (HU) and ||S5} and solv- 
ing for the reduced distance 5" and the image height x' , 
we obtain [8] 

S = . : „ n— = -m x (M n S + M21), (84) 



M12S - M 22 

a; 

~Mi 2 5-M : 



22 



where 



x 

m, = — 



-1 



1 



(85) 



(86) 



a; M 12 ^ - M 22 D 

is the transverse magnification of the system. Notice that 
the reduced output distance S' is a Moebius transform 
of the reduced input distance S. 

4 Image Manifold 

4.1 Height- Angle Phase Space 

a. Object Phase Rectangle. In the height-angle 
phase space, let us construct a rectangle of height dx 



and width d(na): 

fi = xei+naie 2 , (87) 

r 2 = (x + dx)&i + naie 2 , (88) 

fa = (x + dx)ei + (not + d(na))ie 2 , (89) 

r4 = x&i + (na + d(na))ie 2 , (90) 

where ?i, f 2 , fs, and r4 are consecutive vertices of the 
rectangle for a counterclockwise path. (Fig. [6]) 

The two consecutive sides of the rectangle starting 
from f 1 are 

ri2 = h - h = dxei, (91) 

f 2 3 = r 3 - f 2 — d(na)ie 2 . (92) 

The product of these two oriented sides is 

fi 2 f 2 3 = -dxd(na)e 3 , (93) 



because eiie 2 = —63. Notice that Eq. (|93|) is a pure 
vector, whose magnitude is the area of the phase space 
rectangle, as given by Eqs. (|59|) and (f44|) . 




d(na) 



> ei 



Figure 6: An oriented area in phase space denned 
by dxei and d(na)ie2- 



On the other hand, the product of f 12 and r 23 is 

^12^23 = dx d(na)e 3 , (94) 

which is also a pure vector. Notice that the vanishing of 
the scalar component means that the sides of the phase 
rectangle defined by f 12 and f 23 are indeed perpendicular, 
as given by Eqs. gS]) and (|55]) . 

b. Image Phase Rectangle. Since x' and a' are 
functions of x and a, then by the definition of the total 
differential, we have 



dx' 

dx 1 — -^—dx 

ox 



dx 1 



d(n a 



d(na 

d{n'a') d(n'a') 



dx 



-dx 



d{na), (95) 
d{na). (96) 



d(na) 



G 



Using these differential expansions, we can show that the 
image of the phase rectangle defined by the points f± to 
f4 in Eqs. ([87]) to (|90|) is parallelogram defined by the 
following points: 



x ei + n' a' ie^, 



r 4 



(97) 

= jx' + ^ dx)ei + {n'a' + ^ dx)ie 2 , (98) 

dx^ dx' 

— jx' + — — dx + —, r d(na)) ei 

ox djna) 

, , d{n'a>) djn'a') 
+ jn a H + — — r 

OT djna) 

dx' 

= ( x ' + ^7 — \ ^M) e i 

a(na) 



d(na))ze 2 , (99) 



. , . djn a ) ... 
+ (n a + -TT7 — r- a(na)ke 2 . 



(100) 



Notice that dx — at points f\ and r4, and d(na) = 
at points r\ and r-z; both and d(na) are nonzero at 
point f 3 . (Fig. 

The two consecutive sides of the parallelogram image 
starting from ?i are 



dx' , 
— dxex 
ox 



djn'a') 
dx 



dx ie 2 , 



(101) 



' 23 



dx' ,, . d(n'a') . , „„. 
d(na) e 1 + -±- — d{na) ie 2 . (102) 



d{na) v ' d(na) 

Their product is a scalar-vector cliffor: 

'da? dx' d(n'a')d(n'a') 
dx d{na) dx d{na) 

fdx'd{n'a') dx' d(n'a') 
3 \ dx d(na) d(na) dx 



12 ' 23 



dx d(na) — 

dx d(na). 
(103) 



Notice that the magnitude of the vector part is the area 
of the parallelogram in the image phase space. The ratio 
of this area with that of the object in the vector part of 
Eq. IJOSJ is 



dx' djn'a') dx' d{n'a') 
[x ' na] ^ a = -^-d^a)-d(na) dx ' (1 ° 4) 

which is the Poisson commutator bracket. We shall later 
show that this commutator is an invariant [20] that is 
equal to unity. 

On the other hand, the product of f' 12 and r 2 \ is 



12 '23 



e.i 



dx' dx' + djn'a') djn'a') 
dx d(na) dx d(na) 
dx' djn'a') dx' djn'a') 
dx djna) djna) dx 



dx djna) — 

dx djna). 
(105) 



Since the scalar part does not vanish, then the sides f 12 
and r23 of the image parallelogram in phase space are not 
perpendicular. For convenience, let us write the scalar 
part coefficient of dx djna) as 

'dx' dx' djn'a') djn'a') 



, n a }x,na — 



dx djna) 



dx djna) 



(106) 

We shall refer to this as a Poisson anticommutator 
bracket. 



n a 




Figure 7: An oriented area in phase space defined 
by dx' — r' 2 — r[ and idjn'a') = r' 3 — f' 2 . 



c. Symplectic Condition. The differentials of 
Eqs. ([65j) and ((66} are 

dx' = Adx + Bdjna), (107) 

djn'a') = -Cdx + Dna. (108) 

Comparing these equations with Eqs. (|95[) and 
arrive at 



.4 



-C 



dx[ 
dx ' 
djn'a') 
dx 



B = 



D 



dx' 
djna) ' 
djn'a') 

djna) 



(109) 
(110) 



Except for the sign of C, these equations are the known 
partial derivative expressions for the system matrix 
parameters. [21] 

Using the relations in Eqs. (|109p and (|110p . the Poisson 
bracket expressions in Eqs. (|106|1 and (|104|) simplifies to 

{x',n'a%, na = -jAB-CD), (111) 

[x\ria'] x , na = AD + BC=1, (112) 
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where we used Eq. (|62|) . The first equation states that 
in the phase space, image parallelogram does not pre- 
serve the perpendicularity of the sides of the original 
object parallelogram, while the second equation states 
that the object and image parallelograms have the same 
area. This last condition means that the transformation 
defined in Eqs. ([65]) and (|66[) is symplectic. 



4.2 Distance-Height Phase Space 

a. Object Rectangle. Let us define a rectangular 
longitudinal object of length dS and height dx: 

ri = —Sea+xei, (113) 

r 2 = {-S + dS)e 3 + xe 1 , (114) 

r 3 = (-S + dS)e 3 + (x + dx)e u (115) 

r 4 = —Se 3 + (x + dx)ei, (H6) 

where ri, r 2 , r3, and r 4 are the four consecutive vertices 
of the rectangle for a counterclockwise path. (Fig. [5]) 




> e 3 



Figure 8: An oriented area defined by dS e3 and dxei 



The two successive sides of the rectangle are 



ri2 = r 2 - ri = dSe 3 , 
r 2 3 = r 3 -r 2 =da;e 1 . 



Their product is 



r 12 r 23 = dS dxe 3 ei, 



so that 



r 12 • r 23 — 0, 

ri2 Ar 23 = dSdxe 3 ei. 



(117) 
(118) 



(119) 



(120) 
(121) 



Thus, the sides of the rectangle are perpendicular and 
the area of the rectangle is dS dx. 



a. Image Parallelogram. The differential of r' in 
Eq. (|68]) is 

dr' = dS'e 3 + dx'ei, 



where 



, OS OS 
= ~0S~ ~0x 

dx' = ^rrrdS + ~^—dx, 

OS Ox 



(122) 

(123) 
(124) 



by the definition of a total differential. 

Taking the partial derivatives of S' and x' in Eqs. ([84 
and (|85[) with respect to S and x, we get 



OS' 

~os 

0& 
Ox 

~os 

0x[ 
Ox 



(M 12 S - M 22 y 



M 12 x 



{M 12 S - M 22 f 

1 

(M 12 S - M 22 ) " 



-m% (125) 
(126) 

-m 2 x Mi 2 x, (127) 
(128) 



where we used the unit determinant property of M^ ox in 
Eq. (|75[) and the relation for the transverse magnification 
m x in Eq. ([86| . Notice that Eq. (|125[) is the same as the 
expression for longitudinal magnification for infinitesimal 
longitudinal displacement [22j . 

Now, the image of the points ri, r 2 , r 3 , and r 4 in 
Eq. (fTT3|) to (fTT6)l are 



r[ = S'es + x'e-L, 



(129) 



(S' + ^dS)e 3 + (x' + ^dS) ei , (130) 
, n i OS OS . 

{s + Js dS + ^ dx)e3 

0x f Ox' 
+ (x' + —dS + —dx) ei , (131) 

OS' Ox' 
r 4 = (S' + — dx)e 3 + (x' + — dx)e u (132) 
Ox Ox 

respectively. Notice that these points now define a par- 
allelogram. (Fig. [9]) 

The two successive sides of the rectangle are 



L 12 



OS' _ Ox' „ 
—dSe^-dSe,, 



OS' 



Ox' 



1-23 = ~0x dxe3 + ~0x dxei 



Their product is 



r 19 r 9. 3 — r 12 ' r 23 "I" r 12 ^ r 23 



(133) 
(134) 

(135) 
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where 




1 1 > e 3 

S' 

Figure 9: The oriented area defined by dx! = 
f' 2 — f[ and id(n'a') = r' s — f' 2 . 

The coefficient of dS dx in Eqs. (|136p and (|137[) may- 
be expressed in terms of the Poisson anticommutator and 
commutator brackets as 

r „. dS'dS' dx' dx' 

&>^ 8 * = Wte + eslb' (138) 

n dS' dx' dx' dS' 

= ^~dS^ (139) 

Using the magnification equations in Eqs. (|125|) to (|128D . 
Eqs. (fl38]) and (fT39]) becomes 

{S',x'} s , x = -m s x M 12 x, (140) 
[S',x'] s . x = -m 3 x . (141) 

The perpendicularity measure in Eq. (|140p states that the 
image of a longitudinal rectangle will also be a rectangle, 
provided that M\i = 0, which is the characteristic of 
telescopic systems [23], And the area ratio in Eq. (|141[) 
states that the area of the image is proportional to that 
of the object by a factor of negative of the cube of the 
magnification. 



5 Conclusions 

We used geometric algebra to compute the inner (dot) 
and outer (wedge) product of two vectors. We used the 
former to define the Poisson commutator bracket for mea- 
suring the perpendicularity of two vectors; the latter for 
the Poisson anticommutator bracket for measuring areas. 
We adopted the complex height-angle ray formalism we 
developed in previous papers to write down the 2x2 
matrix equations for tracing. And from these equations 
we derive the partial Moebius transforms that relates the 
height and distance of the input ray to that of the output 
ray, as measured from the input and output sides of the 
optical black box. 

For the case of distance- height rays, we define the ob- 
ject to be a differential rectangle in the zx— plane. We 
showed that its image does not preserve the area of the 
object nor the perpendicularity of its sides. A special 
case is that of telescopic systems, where the perpendic- 
ularity of the sides is preserved but not the area of the 
image which is equal to the negative of the cube of the 
transverse magnification. The negative sign means that 
the orientation of the areas is opposite: one area is ori- 
ented clockwise; the other counterclockwise. 

For the case of the height- angle rays, a similar com- 
putation is not possible because we adopted a complex 
height-angle ray formalism, where height is along the 
x— axis and the angle is along the imaginary y— axis. To 
define the pependicularity of two complex vectors, we 
took the scalar part of the product of the complex height- 
angle ray vector with its spatial inverse: if this scalar 
part is zero, then the two complex vectors are perpendic- 
ular. To define the area of the phase space parallelogram 
formed by two complex vectors, we took the magnitude 
of the vector part of the product of the two complex vec- 
tors. (These procedures are only valid if the vector and 
imaginary vector components of each vector are perpen- 
dicular in the complex sense.) We showed that the area 
of differential height-angle rectangular objects in phase 
space is preserved by a 2 x 2 matrix transformation; how- 
ever, the sides of the rectangle are generally not anymore 
perpendicular. The ratio of area of the image in phase 
space to that of the object is equal to the magnitude of 
the Poisson commutator bracket. 

In Differential Geometry, the Poisson commutator 
bracket is already used and well-known. We also hope 
that its anticommutator counterpart would also be sim- 
ilary used, within the context of geometric algebra. 
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